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Abstract 

We study statistical properties of two-dimensional turbulent flows. Three systems are considered: the 
Navier-Stokes equation, surface quasi-geostrophic flow, and a model equation for thermal convection in the 
Earth's mantle. Direct numerical simulations are used to determine 1-point fluctuation properties. Comparative 
study shows universality of probability density functions (PDFs) across different types of flow. Especially for 
the derivatives of the "advected" quantity, the shapes of the PDFs are the same for the three flows, once 
normalized by the average size of fluctuations. Theoretical models for the shape of PDFs are briefly discussed. 
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The central idea of classical turbulence theory is that 
certain statistical properties in turbulent flow are inde- 
pendent of the details of the flow, like its boundaries, 
dissipation mechanism, and the kind of forcing, as long 
as the Reynolds number is sufficiently high. In this sense 
turbulent flow would be universal. Here we shall inves- 
tigate independence not of boundary conditions, forcing 
etc., but look for universality across equations. This 
is not a crazy idea at all and has long been demon- 
strated for other classes of partial differential equations. 
Many partial differential equations with wave-like be- 
havior, possessing common symmetries, exhibit iden- 
tical fluctuation properties, once normalized by their 
standard deviations. (The Schrodinger equation belongs 
to this class, as formulated in the well-known Bohigas- 
Giannoni-Schmit conjecture). In the present paper it 
is demonstrated that three different hyperbolic partial 
differential equations (which all possess the same sym- 
metries) exhibit, under analogous conditions, the same 
statistics for their fluctuations, once normalized by their 
standard deviations. 

The three flows are described by advection-diffusion 
equations 

— +v-Ve = DV^e + f. (la) 

The scalar quantity advected is 6. The forcing / supplies 
the energy dissipated via a dissipation constant D. The 
velocity w is a function of 9, best written in fourier space, 

^ - . k X §(k) 

The two-dimensional cross product fc x is to be under- 
stood as a vector of length \k9\ and direction perpendic- 
ular to k. It follows that V • w = 0. 

Different values of a correspond to different flows [|] . 
The two-dimensional Navier-Stokes equation is a = 2 
and corresponds to the vorticity V x u . The surface 
quasi-geostrophic equation, a — 1, is a special case of 
the important quasi-geostrophic equation that describes 
flow of a shallow layer on a rotating sphere, as relevant 
for planetary atmospheres and oceans ||^. In this case, 
9 is physically interpreted as temperature, which drives 
the flow through its buoyancy effect. The third equation 
considered is a = 3 which also appears in geophysical 
context as a limiting case of a shallow flow on a rotating 
sphere with uniform internal heating Q. Also here, 9 is 
a temperature. 

Other values of a, integer or not, could be consid- 
ered, but this is not done here. Finite-time singularities 



develop for a < 1, but not for a > 1. For the bor- 
der case of a = 1 no corresponding analytical proof has 
been achieved, but numerical evidence suggests that no 
finite-time singularities occur ||] . 

Multiplying eqn. ( pa| ) with 9 and averaging over 
space with periodic boundary conditions yields 

^^{9^) = -D{\V9\^) + {f9). (2) 

Consequently the left hand side of eqn. (|l^) conserves 
(0^) for all a. For a — 2 also (i/^) is a conserved quan- 
tity, while (tJ^) = ^61^) for a = 1, and (-S^) is not con- 
served for a — 3. Equation (1) is invariant under r — )■ 
—r as well as the set of simultaneous transformations 
r rX, t ^ tX^, 9 6'/A", and / ^ //A"+2 (this is 
essentially the Reynolds number invariance) . The family 
of flows with different a has been named a-turbulence 

although this term appears in the literature also for 
other kinds of flow. 

The flow is simulated in a doubly periodic square 
box. Forcing acts on large scales, 4 < |/c| < 6, with 
constant amplitude but random phases renewed at each 
time step. (The time step is constant). Two-dimensional 
turbulent flows produce vortices that merge and grow 
ever larger. These vortices must be destroyed in order 
to reach an equilibrium state. This is done by adding a 
large scale dissipation —76, with 7 small, to the right- 
hand side of eqn. (|l^), restricted to < < 3. All 
these conditions are designed to produce isotropic and 
homogeneous flow. The simulations are carried out with 
a pseudo-spectral method over long time periods using 
4-th order Runge-Kutta integration. A mild spectral 
filter has been used, without complete dealiasing, since 
it is not clear whether complete dealiasing improves or 
worsens the quality of simulations. 

The aforementioned invariance naturally defines a 
Reynolds number for flow of any a as Re = UL/D, 
where U and L are a velocity and length scale respec- 
tively. We choose U — y/ (iP) and L ~ 1 for a large- 
scale Reynolds number. With this definition the max- 
imum Reynolds numbers achieved are on the order of 
several thousands on a 1024x1024 grid for each of the 
three fiows. 

In this paper only 1-point probability density func- 
tions (PDFs) are studied. First, the Navier-Stokes equa- 
tion (a = 2) is treated, which is important by itself and 
also exemplifies the variations and dependencies in the 
PDFs within one equation. In the later part equations 
with different values of a are compared with each other, 
which is the central concern of this paper. 
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The PDFs are obtained from spatial snapshots of the 
flow field and then averaging over 8-24 such ensembles. 
Some of the PDFs shown are scaled by their average 
fluctuation, defined as 



dx\x\P(x). 



(3) 



The integral is over all x. Instead of the first moment the 
standard deviation could have been used equally well. 
All PDFs of the Navier-Stokes equation presented here 
agree with the ones reported from recent simulations by 
Takahashi and Gotoh Q at higher Reynolds numbers. 

Fig, p] shows PDFs for different Reynolds numbers. 
In Fig. llR we see similar but not at all identical shapes 
for the PDFs, a behavior representative for the PDFs of 
other quantities as well. 

According to Fig. ^a velocity components are dis- 
tributed Gaussian. The PDFs for and Vy are almost 
identical, as must be true for isotropic turbulence. The 
tails show the same behavior as seen in decaying two- 
dimensonal turbulence. For a detailed discussion of this 
issue see , where the tails are explained from the influ- 
ence of large vortices, which cause the largest velocity 
gradients. If the two velocity components are statis- 
tically independent of each other, then the PDF of the 
absolute value of v should be a two-dimensional Maxwell 
distribution 



P{x) = ^ exp - 



X 

2^ 

The parameter s is thereby the standard deviation of the 
Gaussian distribution for Vx- The Maxwell distribution 
plotted as dotted curve in Fig. |^b, hence contains no 
free parameter. It is a good first-order approximation. 

As a matter of space not all PDFs can be presented 
here. The scalar (vorticity) is Gaussian in the center. 
The velocity derivatives are also Gaussian. This is par- 
ticularly striking, since velocity derivatives of decaying 
turbulence are not Gaussian |7j. Their core behaves 
much more like a Cauchy distribution 



P{x) 



which has an inflection point even on a logarithmic plot. 
A Cauchy distribution follows theoretically from a "di- 
lute gas" of point vortices of equal strength that move 
randomly, see H. For forced two-dimensional turbu- 
lence dxO (and v^6) could be interpreted as Cauchy dis- 
tribution, but point vortices cannot make any sensible 
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Figure 1: Probability density functions for the Navier- 
Stokes equation at different Reynolds numbers a) vor- 
ticity derivative dxd and b) vorticity dissipation D\V6\'^. 
Many small fluctuations account for most of the dissi- 
pation. The insets show the same data on a logarithmic 
scale. 

predictions on vorticity derivatives. We shall return to 
velocity derivatives later. 

PDFs in a-turbulence have been previously reported 
in (Fig. 7), where a "remarkable similarity" of the 
PDFs for 9 has been pointed out. 

Fig. H shows PDFs for different types of flow. In each 
figure the PDFs of all three fiows are shown simultane- 
ously, and the differernt figures show scalar derivatives 
dxO, scalar dissipation £'|V6'p, and velocity components 
Vx- Apparently the PDFs for the different fiows are the 
same. The agreement is for small as well as large fiuc- 
tuations up to several standard deviations. The very 
largest fiuctuations are inevitably undersampled, which 
accounts for deviations in the tails, which are likely to 
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Figure 2: Probability density functions for vefocities of 
the Navier-Stokes at Re=4500. a) velocity component 
Vx and b) absolute value of velocity \v\. The dotted lines 
are theoretical fits (Gaussian and Maxwellian) . 



fall within measurement errors. 

PDFs by themselves are known not to be particularly 
robust or universal, and it is hence important to compare 
flows under analogous conditions. The Reynolds num- 
bers in the simulations for a = 1,2,3 were Re = 3900, 
4500, 4200 respectively. The differences do not appear 
significant. Comparisons at a lower Reynolds number 
(around 1300) yield the same universalities. 

Also (not shown) overlaps with the same ac- 
curacy as dxd and jVSp. But not all PDFs overlap as 
accurately as the derivatives of 9. The deviations in Fig. 
||c are somewhat larger. Other PDFs show even larger 
deviations, but in none of the PDFs investigated was 
there any drastic difference. These deviations could lie 
within measurement errors. It turns out that quanti- 
ties with larger deviations also show greater variability 
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Figure 3: The main result. Probability density functions 
for different a. The flows have the same forcing and 
similar Reynolds number, a) the gradient dxO b) the 
scalar dissipation D\\76\'^, and c) the velocity component 
Vx- The PDFs are divided by their average fluctuation 
(7, except for b, where this is not necessary. The shape 
of the PDFs is independent of the type of flow. 
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within the same equation for different Reynolds num- 
ber. It appears that some PDFs simply require longer 
averaging than others. It is consistent with the data 
that all 1-point PDFs are identical, but it is also con- 
sistent that some are not. At the very least, PDFs for 
dxd, |V0|^, V^0, Vx (and dyO, Vy) agree within several 
standard deviations. 

We note (from Figures |^a and |3^) that the variation 
in the PDFs for different a is smaller than the variation 
with Reynolds number. It is astonishing that universal- 
ity across equations is more robust than within the same 
equation at different Reynolds numbers. 

The three quations describe physically quite different 
flows and also their mathemtical properties are diverse. 
For example, the velocity is a local function of the scalar 
only for one of the equations, and the kinetic energy is 
conserved in only two of the three flows. Inspite of these 
vast differences, their fluctuation properties differ only 
in absolute size. 

PDFs of velocities on the ocean surface have recently 
been measured |^ . Although none of the equations de- 
scribed here is directly applicable to this situation, the 
velocity components and their derivatives (the latter are 
not shown here) apparently show the same behavior 
as forced a-turbulence. Universality across equations 
might explain the observed PDFs. 

Several, if not all, probability density functions of 
three different flows, described by hyperbolic partial dif- 
ferential equations with reflection symmetry and Reyn- 
olds number invariance, have identical shapes. This 
clearly demonstrates that fluctuation properties can be 
explained in terms of statistical considerations only, and 
do not (or only very marginally) involve the dynamics 
of the flow. 
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